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Rotational evolution of young pulsars due to superfluid de- 
coupling 
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Pulsars are rotating neutron stars that are seen to slow down, and the spin-down rate is thought to 
be due to magnetic dipole radiation. 1,2 This leads to a prediction for the braking index n, which is a 
combination of spin period and its first and second time derivatives. However, all observed values 3 of 
n are below the predicted value of 3. Here we provide a simple model that can explain the rotational 
evolution of young pulsars, including the n — 2.51 of the 958-year-old pulsar in the Crab nebula. 4 The 
model is based on a decrease in effective moment of inertia due to an increase in the fraction of the 
stellar core that becomes superfluid as the star cools via neutrino emission. The results suggest that 
future large radio monitoring campaigns of pulsars will yield measurements of the neutron star mass, 
nuclear equation of state, and superfluid properties. 



The core of a neutron star has densities near and above 
nuclear saturation and extends to ninety percent of the 
radius of the star; the remaining kilometer or so is the 
stellar crust. The core is composed of degenerate mat- 
ter, mostly neutrons and a small fraction of protons and 
electrons (and possibly exotica, such as hyperons and de- 
confined quarks, which we do not consider here). Im- 
mediately after neutron star formation, this matter is in 
a normal state due to the high temperatures reached in 
stellar core collapse. However, neutron stars cool rapidly 
through the emission of neutrinos, and when the temper- 
ature drops below the (density-dependent) critical tem- 
perature for Cooper pairing, neutrons and protons form 
a superfluid and superconductor, respectively. 5,6 Super- 
fluid neutrons rotate by forming quantized vortices, and 
the spatial distribution of these vortices determines the 
rotation rate of the superfluid core, e.g., vortices migrate 
away from the stellar axis of rotation when the superfluid 
angular velocity fi s f decreases while J7 s f cannot change 
if the vortices are fixed in location, i.e., when they are 
pinned. Meanwhile, normal matter (e.g., in the crust) ro- 
tates at an angular velocity ft that decreases as a result of 
energy loss from the stellar surface due to magnetic dipole 
radiation, i.e., dE/dt = -/3ft 4 , where /3 w B 2 R e /6c 3 and 
B and R are the neutron star magnetic field and radius, 
respectively. 1,2 

A rapid decline in surface temperature was recently de- 
tected in the 330-year-old neutron star in the Cassiopeia A 
supernova remnant. 7,8 The observed cooling can be un- 
derstood as being caused by the recent onset of neutron 
superfluidity in the core of the star, combined with a much 
earlier onset of proton superconductivity. 9,8 This has pro- 
vided the first direct constraints on core superfluid and su- 
perconducting properties from neutron star observations. 



These new results motivate studies of possible implica- 
tions. Here we explore the rotational evolution of young 
pulsars using the newly-constrained superfluid properties 
and assuming this superfluid core is allowed to decouple 
(as discussed below). We use simulations 10 of the cooling 
of a neutron star to determine the fraction of the neutron 
star core that is superfluid as a function of time; this al- 
lows us to track the normal and superfluid components of 
the moment of inertia as the star ages (see Figure 3 in 
Supplementary Information) . 

We consider a simple phenomenological model for rota- 
tional evolution of the normal and superfluid components 
of the star: 
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where I and I s { are the moments of inertia of the nor- 
mal and superfluid components, respectively, and A p ; n 
and A m f are torques associated with vortex pinning 11 and 
dissipative mutual friction, 12 respectively. Note that it is 
the rotation of the normal component SI that is observed 
in pulsars. There are three simple limits that we can con- 
sider. The first is that friction acts on a much shorter 
timescale than the spin-down timescale 12 ; this is the con- 
ventional view of rotational evolution, which leads to f2 s f 
closely tracking SI and a braking index n = 3, at odds 
with all measured values. 3 The second limit is when there 
is no pinning or friction; we find that this leads to n > 3. 
The final case, which we consider in detail here, is when 
pinning causes Sl s f ps (x is the first time derivative of 
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parameter x). The evolution equations can then be com- 
bined to give 
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The spin lag, il s f — is the difference in rotational veloc- 
ity between the superfluid and normal components. Some 
examples of decoupled spin evolution are shown in Fig- 
ure 1. 

Conventional pulsar spin evolution only accounts for 
the second term on the right hand side of eq. (3), with 
constant B and I. In this case, pulsars born at a particu- 
lar spin period (P — 27r/f2) evolve by moving along lines 
of constant 13,14 B, and the characteristic age r c (= P/2P) 
is an estimate of the true age of the pulsar. However, this 
again suggests that the conventional picture is incomplete: 
in cases where an independent age can be estimated (e.g., 
from studying the expansion of an associated supernova 
remnant), the result is often quite different from the char- 
acteristic age; this can be seen in Table 1. When su- 
perfluid decoupling is taken into account, spin evolution 
is similar to the conventional course, except now the mo- 
ment of inertia decreases over time. If the spin lag remains 
small (e.g., as a result of an angular momentum sink act- 
ing on the superfluid), only a small deviation (from evo- 
lution along a constant B track) is seen at intermediate 
and late times (after ~ 1000 — 2000 yr). Note that we are 
considering intermediate times (<~ 10 3 — 10 5 yr) in the life 
of a neutron star, in between the short timescale for glitch 
recurrence 15 (~ 1 yr) and the long timescales for magnetic 
field diffusion 13,14 and cooling 16 - 18 (~ 10 5 - 10 6 yr). On 
the other hand, if the spin lag is allowed to become large 
(e.g., due to strong vortex pinning), then we see that a de- 
creasing / can mimic a strongly increasing magnetic field. 
A further departure from the conventional picture is that 
the characteristic age is not an accurate indication of the 
true age of a pulsar at early times. 

Not only are the spin period P and first time derivative 
of the period P observable quantities, in some cases the 
second time derivative P can be measured (x is the second 
time derivative of parameter x). This provides another 
test for our model. Table 1 gives data on eight systems 
where the braking index n, which is proportional to P, 
has been measured, while our model predicts the braking 
index to be [see eq. (3)] 
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When f2 s f — fl <C O, the last term is irrelevant, and the 
braking index is simply given by n = 3 — 4t c \I/I\. 

Since r c and n are observable quantities (related sim- 
ply to P, P, and P) for a given pulsar, we can compare 
our predictions to the pulsars in Table 1; this is shown in 
Figure 2. For the Crab pulsar (the only one with a known 
age), we infer a mass of ss 1.8Ms un (for the particular 



equation of state and superfluid properties we consider; 
see Supplementary Information). Furthermore, we can 
use the mass determination to estimate the initial period 
and magnetic field of the pulsar, and we find an initial 
period - 0.02 s and B~4x 10 12 G. 

The fact that our simple model is able to explain the 
observed pulsar properties demonstrates the merits of the 
notion of decoupled spin-evolution due to the onset of core 
superfluidity. However, key questions remain to be un- 
derstood. The main assumption in our model is that core 
superfluid neutrons are allowed to decouple and pin. Core 
pinning is thought to be the result of the interaction be- 
tween superfluid vortices and fluxtubes in the proton su- 
perconductor. 11 Whether this mechanism can be strong 
enough to act in the way assumed in our analysis is not 
clear at this time. Theoretical work is also required to 
determine whether the spin-lag O s f — f2 can be kept small 
during the evolution by an (at this time) unspecified an- 
gular momentum "sink" . From an observational point-of- 
view, the model is promising, in particular since it can be 
used to infer the mass of individual pulsars. In the fu- 
ture, discovery and long-term monitoring of a large num- 
ber of systems by radio telescopes such as LOFAR and 
SKA will allow accurate timing of many young pulsars. 
Taking these pulsars as an ensemble, we can constrain 
the nuclear equation of state and superfluid properties 
since these determine the evolution of the moment of in- 
ertia. Finally, radio timing measurements may be able to 
constrain neutron star thermal evolution, independent of 
measurements at X-ray energies. 

Received 17 August 2012; accepted — . 

1. Pacini, F. Rotating neutron stars, pulsars and supernova 
remnants. Nature 219, 145-146 (1968). 

2. Gunn, J. E. & Ostriker, J. P. Magnetic dipole radiation 
from pulsars. Nature 221, 454-456 (1969). 

3. Espinoza, C. M., Lyne, A. G., Kramer, M., Manch- 
ester, R. N., & Kaspi, V. M. The braking index of 
PSR J1734— 3333 and the magnetar population. Astro- 
phys. J. 741, L13-1-6 (2011). 

4. Lyne, A. G., Pritchard, R. S., & Graham-Smith, F. 23 
years of Crab pulsar rotational history. Mon. Not. R. As- 
tron. Soc. 265, 1003-1012 (1993). 

5. Migdal, A. B. Superfluidity and the moments of inertia of 
nuclei. Sou. Phys. JETP 10, 176-185 (1960). 

6. Baym, G., Pethick, C, & Pines, D. Superfluidity in neu- 
tron stars. Nature 224, 673-674 (1969). 

7. Heinke, C. O. & Ho, W. C. G. Direct observation of the 
cooling of the Cassiopeia A neutron star. Astrophys. J. 
719, L167-L171 (2010). 

8. Shternin, P. S., Yakovlcv, D. G., Heinke, C. O., Ho, 
W. C. G. k, Patnaude, D. J. Cooling neutron star in the 
Cassiopeia A supernova remnant: Evidence for superfluid- 
ity in the core. Mon. Not. R. Astron. Soc. 412, L108-L112 
(2011). 



Young pulsars and superfluidity 3 



Pulsar Supernova Period Period derivative Characteristic Braking Age 

Name remnant (s) (s s _1 ) age r c (yr) index n (yr) 



B0531+21 


Crab 


0.0331 


4.23x10" 


13 


1240 


2.51(1) [4] 
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J0537-6910 


N157B 


0.0161 


5.18x10" 


14 
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-1.5(1) [20] 
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B0833-45 
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1.4(2) [24] 
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[25] 


J1119-6127 


G292.2-0.5 


0.408 
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1610 


2.684(2) [26] 


nool^c, 


[27] 


B1509-58 


G320.4-1.2 


0.151 


1.54x10" 


12 


1550 


2.839(3) [22] 


< 21000 


28] 


J1846-0258 


Kesteven 75 


0.325 


7.08x10" 


12 


729 


2.65(1) [22] 


10001?°°° 


[29] 


J1734-3333 


G354.8-0.8 


1.17 


2.28x10" 


12 


8120 


0.9(2) [3] 


> 1300 



Table 1. Pulsars with observed braking index. Periods and period derivatives are taken from ref. 19. Numbers 
in parentheses show braking index uncertainty in the last digit. For J1734— 3333, we give a lower limit of the age, 
which we estimate by (1) considering the supernova remnant size (21 parsecs; ref. 30) and remnant expansion velocity 
i'Snr, in order to obtain an age ~ 2000 yr (10 4 km s~ 1 /«snr) and (2) considering the pulsar's distance away from 
the centre of the supernova remnant (46 parsecs; ref. 30) and pulsar space velocity i> pu isar, in order to obtain an age 
~ 23000 yr (2000 km B-Vwpuhar). 
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Figure 1. Pulsar spin period versus spin period 
derivative. Open circles are observed values taken from 
the ATNF Pulsar Catalogue. 19 Stars denote pulsars with 
a measured braking index (see Table 1). Thin diag- 
onal lines denote characteristic age (= P/2P) and in- 
ferred magnetic field [= 3.2 x 10 19 G {PP) 1 ' 2 ]. Curves 
are spin evolution tracks for pulsars with mass 1.8Ms un 
(red, dashed) and 1.4Mg un (blue, solid), where the spin 
lag is maintained at O s f/fi — 1 < 10 -6 (by an additional 
angular momentum sink); filled circles denote the evolu- 
tion at ages 1000, 2000, and 10 4 yr. Also shown (green, 
thin solid) is the evolution of a pulsar where the spin lag 
is allowed to grow. From left to right, the initial spin 
period and magnetic field (P,B) are taken to be (0.02 s, 
5 x 10 12 G), (0.1 s, 10 13 G), and (1 s, 10 14 G). 




Age (yr) 

Figure 2. Constraining neutron star properties. 

Crosses denote pulsars with a measured braking index and 
an estimate or limit on their true age (see Table 1). Curves 
show the evolution of the effective moment of inertia 
for pulsars with mass 1.8Ms un (short-dashed), 1.4Ms un 
(solid), and 1.2Ms U n (long-dashed), assuming a particu- 
lar stellar and superfluid model, 9 while the shaded region 
represents an alternative model 8 ; both sets of superfluid 
parameters fit the rapid cooling seen in the Cassiopeia A 
neutron star (see Supplementary Information). These mo- 
ment of inertia evolution curves are analogous to the ther- 
mal cooling curves 16 " 18 used to determine properties of 
neutron stars. 
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Supplementary Information 

Neutron star cooling and moment of inertia evo- 
lution 

Neutron stars begin their lives very hot (with tempera- 
tures T > 10 11 K) but cool rapidly through the emis- 
sion of neutrinos. Neutrino emission processes depend 
on uncertain physics at the supra-nuclear densities (p > 
2.8 x 10 14 g cm -3 ) of the neutron star core. 16-18 The re- 
cent observation of rapid cooling 7 ' 8 of the neutron star 
in the Cassiopeia A supernova remnant provides the first 
constraints on the (density-dependent) critical tempera- 
tures for the onset of superfluidity of core protons T cp (in 
the singlet state) and neutrons T cnt (in the triplet state), 
i.e., T cp ~ (2-3) x 10 9 K and a maximum T cnt w 5 x 10 8 K 
for a superfluid neutron model at relatively shallow densi- 
ties 9 or a maximum T cnt « (7 — 9) x 10 8 K for a superfluid 
neutron model at deep densities. 8 

We calculate the thermal evolution of neutron stars 
by solving the relativistic equations of energy balance and 
heat flux. We use a stellar model based on the APR equa- 
tion of state and consider either the shallow or deep model 
for triplet neutron pairing in the core. 10 Note that we 
improve upon the calculations of ref. 10 by using more 
accurate ion and electron heat capacities. 31 Due to high 
thermal conductivity, neutron stars after ~ 10 — 100 yr 
are essentially isothermal 32-34 ; the exact time is unim- 
portant for our work since the youngest neutron star we 
consider has an age « 1000 yr (see Table 1). The evolution 
of the (gravitationally redshifted) surface temperature for 
models with the shallow neutron superfluid and an un- 
magnetized iron envelope is shown in Figure 3. Enhanced 
neutrino emission due to Cooper pairing of neutrons be- 
gins to occur at the onset of superfluidity in the core. At 
an age of a few hundred to one thousand years, large por- 
tions of the (predominantly neutron) stellar core become 
superfluid, and rapid cooling occurs. 

We calculate the moment of inertia of the normal and 
superfluid components by 35 

1= ^ £( P + P/c 2 )ArUr, (5) 

where P is the local pressure, A = (1 — 2Gm/c 2 r)~ 1 , and 
to is the mass enclosed within r. Note that we neglect 
for simplicity a term that approximately accounts for rel- 
ativistic frame-dragging, 35 which would only change I by 
< 10%. Figure 3 shows the evolution of the normal com- 
ponent of I. 

In order to facilitate numerical calculations of I(t) and 
I(t), we fit I(t) with the following function 

log/ = a 3 tanh[a 2 (logt-ai)] + (logi-6i) b2 +ci, (6) 

where the fit parameters 01,02,03,61,62,01 are given in 
Table 2. A comparison of the approximate fit to the exact 
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Figure 3. Thermal and moment of inertia evo- 
lution of neutron stars. Curves show the evolution 
of the surface temperature (top) and moment of iner- 
tia (bottom) for neutron stars with mass 1.8 Msw ( re d, 
short-dashed), 1.4Ms un (blue, solid), and 1.2Msun (pur- 
ple, long-dashed). Dotted curves show analytic fits using 
eq. (6). 

Neutron star mass ai £12 03 61 62 ci 
1.2 Msun 2790 7A T5 OTTO To 45.24 
1.4 M Sun 2.59 4.6 -1.2 0.95 -1.2 44.97 
1.8 Msun 1.98 2.6 -0.31 -1.6 44.945 

Table 2. Fit parameters for moment of inertia evo- 
lution. 

result is shown in Figure 3. Note that the largest devia- 
tions occur briefly at (early) times that are not of current 
relevance. 



